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ABSTRACT 

We report first results from the "Streaming Motions of Abell Clusters" (SMAC) project, an all- 
sky Fundamental Plane survey of 699 early- type galaxies in 56 clusters between ~ 3000 km s _1 and 
~ 14000 km s _1 . For this sample, with a median distance of ~ 8000 km s -1 , we find a bulk flow of 
amplitude 630±200km s _1 , towards I — 260 ± 15°, b = — 1 ± 12°, with respect to the Cosmic Microwave 
Background. The flow is robust against the effects of individual clusters and data subsets, the choice 
of Galactic extinction maps, Malmquist bias and stellar population effects. The direction of the SMAC 
flow is ~90° away from the flow found by Lauer & Postman, but is in good agreement with the gravity 
dipole predicted from the distribution of X-ray-luminous clusters. Our detection of a high-amplitude 
coherent flow on such a large scale argues for excess mass density fluctuation power at wavelengths 
A ;> 60 /i _1 Mpc, relative to the predictions of currently-popular cosmological models. 

Subject headings: galaxies: distances and redshifts — galaxies: elliptical and lenticular, cD — galaxies: 
clusters: general — cosmology: observations — large-scale structure of Universe 
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of Dale et al. (1998), none of these studies approach the 
depth and sky coverage of the LP survey. 

In this Letter, we report the detection of a significant 
coherent flow from the SMAC survey of cluster peculiar 
motions. The survey comprises 699 early- type galaxies in 
56 clusters mainly within 12000 km s" 1 . A more detailed 
description and analysis of the survey will be presented in 
a forthcoming series of papers. 

2. DATA & METHOD 

The distance indicator used in this study is the Funda- 
mental Plane (FP) of early-type galaxies (Davis & Djor- 
govski 1987; Dressier et al. 1987). The FP relates the ef- 
fective (half-light) radius, logi? e (the distance-dependent 
quantity), the mean surface brightness within this radius, 
(//) , and the central velocity dispersion, logo - . 

The SMAC cluster sample consists of new data for 40 
clusters (to be reported in future papers), supplemented 
with data from the literature for 16 clusters previously 
studied in H97. All but 8 of these are Abell/ ACO clus- 
ters, and all but 5 have czq < 12000 km s _1 . The median 
distance of the SMAC sample is ~ 8000 km s _1 . The num- 
ber of early-type galaxies observed, per cluster, is in the 
range 4 - 56, with a median of 8. 

New velocity dispersions and Mg2 linestrengths were ob- 
tained from the Isaac Newton 2.5m and Anglo- Australian 
3.9m telescopes. We followed the homogenization proce- 
dure of Smith et al. (1997) to bring new spectroscopic data 
and existing data from the literature onto a common sys- 
tem. Corrections of up to 0.03 dex in a are derived for 
each system, but typical uncertainties in these corrections 
are only 0.005-0.008 dex, corresponding to 1.6-2.6% sys- 
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1. INTRODUCTION 

The dipole anisotropy of the Cosmic Microwave Back- 
ground (CMB) radiation is generally interpreted as a 
Doppler effect due to the motion of the Sun with respect 
to the CMB rest-frame. The velocity of the Local Group 
(LG), 627±22kms- 1 towards Z = 276±3°, b = 30 ±3°, 
is now well determined from COBE (Kogut et al. 1993). 
Less well known, however, is the depth and degree to which 
nearby galaxies share in this motion and, by implication, 
the scales and amplitudes of the mass fluctuations respon- 
sible for the flow. 

To a depth of ~ 6000 km s _1 , a rough consensus has 
emerged from recent peculiar velocity surveys of galax- 
ies. For instance, Giovanelli et al. (1998) find a flow 
(in the CMB frame) of 200 ± 65 km s" 1 towards I = 295°, 
b = 25°, from an I-band Tully-Fisher survey, while the 
Mark III velocity compilation yields 370 ± 110 km s _1 to- 
wards I = 306°, b = 13° (Dekel et al. 1998, in prep.). Be- 
yond this depth, however, the situation is much less clear. 
Lauer & Postman (1994, hereafter LP), using a photo- 
metric distance indicator for brightest cluster galaxies, 
found a 689 ± 178 km s _1 bulk motion (towards I = 343°, 
b = 52°) for an all-sky sample of 119 Abell/ACO clusters 
to 15000 km s _1 depth. Currently-popular cosmological 
models have too little large-scale power to generate coher- 
ent flows on such large scales (Feldman & Watkins 1994; 
Strauss et al. 1995; Jaffe & Kaiser 1995). The LP result 
has been challenged by subsequent peculiar velocity sur- 
veys (Riess et al. 1995; Giovanelli et al. 1996; Hudson et 
al. 1997, hereafter H97; Giovanelli et al. 1997; Muller et 
al. 1998; Dale et al. 1998). However, with the exception 
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tematic uncertainty in distance, per observing run. Of the 
699 galaxies in the final SMAC sample, velocity dispersion 
data is drawn wholly from new observations for 41% of the 
sample. For 21% of the sample, the final a derives from 
both new data and published measurements. For the re- 
maining 38%, the dispersions are derived from previously 
published data. 

New R-band photometric data were obtained from the 
Jacobus Kapetyn 1.0m and Cerro Tololo Inter- American 
Observatory 0.9m telescopes. Effective radii and surface 
brightnesses were determined by fitting a one-dimensional 
J? 1 ' 4 -law profile to the circular aperture photometry, with 
corrections for seeing, cosmological effects and Galactic ex- 
tinction (using the maps of Schlegel, Finkbeiner and Davis 
1998; hereafter SFD). Additional photometric data has 
been drawn from the literature and carefully combined 
with newly obtained parameters. Our new photometry 
provides data for roughly half of the 699 galaxies in the 
final SMAC sample. 

Our method for determining cluster distances and bulk 
flows follows H97. We summarize the important points 
here. To obtain cluster distances, we use the inverse form 
of the FP, i.e. we regress on the distance-independent 
quantity loger. This regression is insensitive to photomet- 
ric selection effects (e.g. selection on magnitude or diame- 
ter; see H97, Strauss & Willick 1995). The FP slopes and 
scatter for the SMAC sample are consistent with previous 
results for the inverse FP (H97). The resulting distance 
error is 21% per galaxy, and 3% - 11% per cluster, with 
a median of 7%. We make a correction for Malmquist 
bias under the assumption that clusters are drawn from 
a homogeneous underlying density field. Due to the large 
number of objects per cluster, the Malmquist bias correc- 
tions are small (typically ~ 2%) for this sample. Finally, 
we adjust the FP zero-point so that the cluster peculiar 
velocities have no net radial inflow or outflow with respect 
to the CMB frame. 

A bulk flow model provides the simplest parameteriza- 
tion of the peculiar velocity field. We minimize 

x a = E (T *~r* <)2 a) 

where Vi is the CMB frame peculiar velocity of cluster 
i with direction vector f.;, and V is the bulk flow. The 
total error, ej, is the quadrature sum of the cluster dis- 
tance error (in the range 100-1200km s _1 ), the error in 
the mean cluster cz (typically 150 km s _1 ) and a "thermal 
component" (set to 250 km s^ 1 ) which allows for small 
scale fluctuations around the mean bulk flow. Given the 
error- weighting in equation (Q), the weighted mean depth 
of the sample is ~ 6700 km s _1 . 

We have tested the bulk flow recovery with Monte Carlo 
simulations in which galaxies are assigned peculiar veloci- 
ties consistent with an input bulk flow plus random errors. 
Because our survey has good sky-coverage, we find that 
our bulk flow fit is not affected by "geometry bias", i.e. 
there is little covariance between the monopole and dipole 
(bulk flow) terms. 

3. RESULTS 

The principal result of this Letter is that the SMAC 
sample exhibits a CMB-frame bulk flow of V = 630km s _1 



(error-bias corrected, see LP), towards Z = 260±15°, 
b = — 1 ± 12°. The peculiar velocities of the sample clus- 
ters are shown, as a function of angle from the apex di- 
rection, in Fig. [l[ The linear trend seen here is a clear 
signature of a bulk streaming motion. 

The random error in the bulk flow due to distance and 
velocity errors is 180 km s" 1 . A further uncertainty in the 
bulk flow arises from uncertainties in the corrections to a 
for different runs. We estimate this uncertainty by gener- 
ating bootstrap realizations of the corrections and analyz- 
ing the bootstrap-corrected datasets in the same way as 
the real data. Through this procedure, we find that this 
source of systematic error contributes an uncertainty of 90 
km s _1 to the measured bulk flow. 

The total velocity error in the direction of the flow is 
200 km s _1 . The error ellipsoid is triaxial: the direction 
of smallest error (±100 km s" 1 ) is along I = 314°, b = 49°, 
and the direction of largest error (±214 km s _1 ) is along 
I = 237°, b = -11°, which is within 24° of our bulk flow. 
Allowing for the three degrees of freedom in the bulk flow 
vector we find that our sample is inconsistent with being 
at rest in the CMB frame at the 99.9% confidence level. A 
95% lower limit on the bulk motion is 400 km s _1 . If we 
divide the sample by distance into two parts with equal er- 
rors, the outer shell of the SMAC sample gives a larger am- 
plitude flow (by 240±360km s _1 ) but is consistent, within 
the errors, with that found from the inner sphere. Both 
inner and outer samples independently yield a bulk flow 
with a significance ;> 98%. 



2000 



> 



-2000 



1 


, I I 










1 1 1 1 

A1736 


1 1 














• 


• 


• 






• * 




* • 




< 




• 






. * 

•• ' - . 
















• 




















• 


< 


» 


S030 


i 






•• 


• 














1 


, , 1 













A1139 




-1 


-0.5 









0.5 




1 



COS SMAC 



Fig. 1. — Peculiar velocities as a function of angle from the apex 
of the bulk flow. Symbol sizes are inversely proportional to errors. 
The solid line is the best-fitting bulk flow model, and clusters which 
deviate from this flow at the 2<r level are indicated, with error bars. 

We have explored a wide range of potential systematic 
effects. By excluding each cluster from the sample in turn, 
we find that no individual cluster is responsible for more 
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than 50 km s _1 of the total bulk motion. No single su- 
percluster structure dominates the flow. We find a similar 
result excluding individual spectroscopic runs. Further- 
more, because of the balanced sky coverage, the bulk flow 
is insensitive to the adopted zero-point of the FP relation. 
A 1% change in the zero-point (corresponding to its la 
uncertainty) alters the bulk flow by only ~ 25 km s^ 1 . 
If, instead of SFD extinction corrections, we use those of 
Burstein & Heiles (1982), the recovered bulk flow drops by 
~ 100 km s _1 . In order to assess the effect of Malmquist 
bias and redshift cuts, we have performed a simultaneous 
FP and bulk flow fit ( "Method II" in the terminology of 
Strauss & Willick 1995); we recover a bulk flow which is 
identical within the errors (~ 40km s _1 change). The ef- 
fect of a possible FP dependence on morphological type 
is negligible: E and SO subsamples of the data give con- 
sistent results (they differ by 140 ± 330km s _1 ). Finally, 
we have examined the possible effect of cluster-to-cluster 
stellar population differences, using the Mg2 index as an 
age/metallicity indicator. We find that our clusters are 
consistent with following a universal Mg2- a relation. Fur- 
thermore, we find no correlation between peculiar velocity 
and the (non-significant) offsets of clusters from the Mg2- 
(7 relation. 
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Fig. 2. — Peculiar velocities projected onto a plane between 
SGX=0 and SGZ=0. The negative half of the horizontal axis points 
towards I = 272°, b = —3°, which is only ~ 13° from the direction 
of the SMAC bulk flow. We plot only the 43 clusters within 45° 
of the plane. The shaded region indicates the Zone of Avoidance 
(|6| < 15). The distance to each cluster is indicated by the circle 
and the CMB frame redshift is indicated by the tip of the vector. 
Outflowing clusters have filled circles and solid lines; inflowing ones 
have open circles and dotted lines. Circle size is inversely propor- 
tional to the error in the peculiar velocity. Important locations are 
indicated: the Local Group (LG); the Great Attractor (GA); the 
Shapley Concentration (SC) and the Horologium-Reticulum (HR) 
supercluster. 



4. DISCUSSION 

The motion of the Local Group (LG) and the bulk mo- 
tion of the SMAC sample both lie in the plane which is 
at a 45°angle from the SGX=0 and SGZ=0 Supergalactic 
planes. The CMB-frame peculiar velocities of SMAC clus- 
ters within ±45° of this plane are shown in Fig. ^. This 
figure shows the continued outflow beyond the Great At- 
tractor ("GA", Lyndcn-Bcll ct al. 1988) both above and 
below the Galactic Plane, indicating that the local motions 
are not generated wholly by the GA. On the opposite side 
of the sky the data, though sparser, exhibit an inflow of 
similar magnitude, suggesting a bulk flow with very large 
coherence length. 

The direction of the SMAC bulk flow in Galactic coor- 
dinates is shown in Fig. ||. Also shown are the positions 
of two most prominent concentrations of clusters within 
300 /i _1 Mpc (Tully et al. 1992), namely the Shapley Con- 
centration (Scaramella ct al. 1989; Raychaudhury 1989) 
and the Horologium-Reticulum supercluster (Lucey et al. 
1983). The direction of the SMAC bulk motion is roughly 
between these two concentrations. Indeed, there is some 
evidence in our sample of inflow towards the Shapley Con- 
centration. The data are too sparse in the foreground 
of Horologium-Reticulum to measure any infall. The di- 
rection of our flow agrees very well (< 16°) with the di- 
rection of LG motion predicted from the density field of 
X-ray/Abell clusters (Plionis & Kolokotronis 1998). The 
IRAS PSCz galaxy redshift survey is of sufficient depth to 
yield accurate predicted peculiar velocities (Branchini ct 
al. 1998) for our clusters. A preliminary comparison indi- 
cates excellent directional agreement between the PSCz- 
predicted and the observed bulk flow of SMAC clusters. 

How does our bulk flow compare to results from other 
surveys? One might naively estimate the level of agree- 
ment between two surveys from a \ 2 statistic derived from 
the difference between the two bulk flow vectors and the 
sum of their observational-error covariance matrices. This 
procedure is incorrect, however, because it does not allow 
for the fact that different surveys probe different volumes 
and regions on the sky. As stressed by Watkins & Feldman 
(1995), the incomplete cancelation of small-scale flows in- 
ternal to a volume can have a significant effect on the total 
measured bulk flow. The significance of disagreements be- 
tween flow vectors is always overestimated if these effects 
are neglected. 

For example, according to the naive comparison de- 
scribed above, our result is in conflict with the bulk flow 
vector of LP and of Dale et al. (1998) . Our flow is of similar 
amplitude to that found by LP, but the direction is ~ 90° 
from the apex of the LP dipole, so at face value the agree- 
ment is poor. Dale et al. find a bulk flow of <J 200 km s^ 1 
for their deep Tully-Fisher cluster sample, in contrast to 
our result 5 . A detailed analysis, allowing for the effects of 
survey geometries and internal flows, is required in order 
to determine whether the apparent disagreements with LP 
and Dale et al. are in fact significant. Such an analysis is 
in progress. 

In contrast to the two surveys just discussed, our flow 
agrees very well in amplitude and direction with that of the 



5 It is worth noting, however, that our error-weighted bulk flow is within the 2<r error ellipse of Dale et al.'s volume- weighted solution, 
although not of their error-weighted solution. 
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Fig. 3. — Peculiar velocities and bulk flow directions in Galactic coordinates (Aitoff projection). Clusters in the SMAC survey are indicated 
by open circles (inflowing) or asterisks (outflowing), with symbol size proportional to peculiar velocity. The direction of the SMAC bulk 
flow is indicated by the large "S", while the dots show 1000 directions drawn from the bulk flow error ellipsoid. Important directions are 
indicated: the motion of the Local Group (LG); the motion of the Lynden-Bell et al. (1988) sample (7S); the bulk flow of LP; the motion of 
the Riess et al. (1997) SNIa sample (R); the Shapley Concentration (SC); the Horologium-Reticulum supercluster (HR); and the predicted 
motion of the LG from X-ray clusters (X) (Plionis & Kolokotronis 1998). The bulk flow of Dale et al. (1998) is small and consistent with 
zero. Consequently, its direction is not meaningful and has not been plotted here. 



similarly deep (9000 - 13000 km s" 1 ) Tully-Fisher sample 
of Willick (1998, also private communication). The SMAC 
bulk flow is not inconsistent with preliminary results from 
the EFAR project (Saglia et al. 1998). Finally, we have 
performed a bulk flow fit to the SNIa velocity data of Riess 
et al. (1997), which yields a flow of 320±160km s^ 1 (error- 
bias corrected) directed towards I = 282°, b — —8°. This 
direction is only ~ 20° from the SMAC flow apex. 

Our result can be compared to expectations based on 
popular families of cosmological models. The bulk flow 
amplitude measured by a typical (i.e. randomly-placed) 
observer depends on the power spectrum of mass den- 
sity fluctuations and the geometry of the survey. Follow- 
ing Kaiser (1988), we have calculated the fc-space window 
function for the bulk flow given the SMAC survey geom- 
etry. We find that the SMAC bulk flow probes the mass 
power spectrum on scales larger than A ~ 60/i -1 Mpc. 
If we fix the power spectrum on the largest scales using 
COBE, then the high amplitude of the SMAC bulk flow 
implies substantial power on intermediate scales. For vari- 
ants of the cold dark matter (CDM) cosmology with either 
a cosmological constant, neutrinos, or a tilt (n ^ 1) in the 
initial primordial spectrum, the choice of free parameters 
fixes the power spectrum on sub-COBE scales. For the 
three COBE-normalizcd CDM variants, we require choices 
of the parameters such that the rms fluctuation in mass 
density contrast in an 8 h~ 1 Mpc sphere is erg f2 06 ;> 0.87, 
0.80 or 0.85, respectively, at the 90% confidence level. This 
is in conflict with the value crgfi ' 5 ~ 0.55 inferred from 
the abundance of rich clusters (Vianna & Liddle 1996; Eke 
et al. 1996; Pen 1998), but is in better agreement with the 
value er^ ' 6 ~ 0.8 inferred from other peculiar velocity 
surveys (Kolatt & Dekel 1997; Zaroubi et al. 1997). Alter- 
natively, if the power spectrum is to fit both the abundance 



of rich clusters and the COBE fluctuations, then it must 
be considerably more "peaked" between these scales, rel- 
ative to the CDM variants. 

5. SUMMARY 

We have recently completed a FP survey of 699 early- 
type galaxies in 56 clusters within ~ 12000 km s _1 . For 
this sample, we find a large-scale bulk flow of amplitude 
630 ± 200 km s" 1 towards I = 260 ± 15°, b = -1 ± 12°. 
Our result is robust against the effects of individual clus- 
ters and data subsets, the choice of Galactic extinction 
maps, Malmquist bias and stellar population effects. 

This result suggests that the mass fluctuation power 
spectrum has a high amplitude on scales A ~ 60 — 
600ft.~ 1 Mpc. This regime, near the peak of the power 
spectrum, is at present poorly constrained by other meth- 
ods. On these scales, CMB anisotropies are dominated by 
acoustic oscillations and depend mainly on the baryon con- 
tent, whereas current galaxy redshift surveys are insuffi- 
ciently deep, and are subject to the (unknown) relationship 
of galaxies to the underlying mass density fluctuations. 

Further analyses, currently underway, will (i) determine 
the consistency of our results with those of other surveys, 
allowing for the different sample geometries; (ii) compare 
measured velocities with predictions of the IRAS PSCz 
redshift survey; and (iii) constrain directly the mass power 
spectrum, by means of the dipole and higher moments of 
the velocity field. 

We thank the staff of the Isaac Newton Group of 
telescopes at the Observatorio del Roque de los Mucha- 
chos, the Anglo- Australian Observatory and the Cerro 
Tololo Inter- American Observatory, where these observa- 
tions were conducted. 
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